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Recently, manipulating heat transport by phononic devices has received significant attention, in which 
phonon - a heat pulse through lattice, is used to carry energy. In addition to heat control, the thermal 
devices might also have broad applications in the renewable energy engineering, such as thermoelectric 
energy harvesting. Elementary phononic devices such as diode, transistor and logic devices have been 
theoretically proposed. In this work, we experimentally create a macroscopic scale thermal rectifier based on 
reduced graphene oxide. Obvious thermal rectification ratio up to 1.21 under 12 K temperature bias has 
been observed. Moreover, this ratio can be enhanced further by increasing the asymmetric ratio. 
Collectively, our results raise the exciting prospect that the realization of macroscopic phononic device with 
large-area graphene based materials is technologically feasible, which may open up important applications 
in thermal circuits and thermal management. 

With the development of solid-state devices, the integrated circuits can be built up by diodes and 
transistors. The conventional diodes and transistors can realize "on" and "off state by controlling 
the electrical conduction. If the thermal conduction can be well controlled, we can also realize thermal 
circuits by using these novel phononics devices. In addition to information processing, the thermal devices might 
also have broad applications for heat control/management in the future. Recent years, some theoretical proposals 
for thermal rectifier have been put forward, include structures based on one-dimensional nanotubes and nano- 
wires^"^, and two-dimensional graphenes^^"^^ The first solid-state thermal rectifier has been reported experi- 
mentally by C. W. Chang et aV'^ in carbon nanotubes. 

However, it is still a challenge to create macroscale phononic devices based on these nanoscale architectures. 
Inspired by the study of thermal rectifier in nanoscale, an oxide thermal diode made of two cobalt oxides with 
different thermal conductivities have been experimentally demonstrated^^'^^. The working mechanism for this 
two- segment diode is that they have thermal conductivities with different temperature dependence. From a 
simple classical heat conduction analysis, the rectification of this two-segment diode is proportional to the 
difference in thermal conductivity, and inversely proportional to the interface resistance^^. In most bulk materials, 
the large interface contact resistance between the two segments limits the enhancement of the rectification ratio. 
Thus it is indispensable to develop single component thermal rectifier. 

In this work, for the first time, we demonstrated macroscopic scale thermal rectification phenomenon in 
asymmetric reduced graphene oxide (rGO) materials. According to recent efforts, rCO can be used to fabricate 
paper-like films with large area, excellent mechanical property, high thermal conductivity and electrical conduc- 
tivity^^"^°. It has been demonstrated theoretically that compared to single-layer graphene, a larger rectification 
ratio can be achieved in multi-layer structures, due to the presence of layer-layer interactions^^. Combined with 
the availability of high quality, here we use a "rGO paper" material composed of stacked rGO to make thermal 
rectifier. Benefit from their high thermal conductivity and strong temperature dependence around room tem- 
perature, high thermal rectification ratio is achieved in both triangular and two -rectangular shaped rGO based 
thermal rectifier, which is in good agreement with the numerical calculation. 

Results 

The optical photograph of the rGO paper is shown in Fig. 1(a). Figure 1(b) illustrates the surface of rGO film. 
Some ripples of graphene can be seen clearly. Figure 1(c) shows the cross-section of rGO. The film thickness is 
—50 jim. Figure 1(d) presents the Energy Dispersive X-ray (EDX) spectrum of rGO film. The film is mainly 
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Figure 1 | (a) Photograph of the rGO paper, (b) Surface profile of rGO 
paper under SEM. (c) Cross section image of rGO paper under SEM. The 
thickness is —50 um. (d) EDX spectrum of rGO paper. There remains fittle 
oxygen atoms, (e) Schematic view of the measure system. A triangular rGO 
structure is tested, (f) Schematic view of the devices. 

graphene stacked layers with little oxygen remained. The ratio 
of carbon and oxygen is —5.56. The precise chemical composi- 
tion and structure of rGO film can be found as the Supplementary 
Fig. SI —S3 online. Our measurement system is described as a schem- 
atic figure in Fig. 1(e). The rGO samples are suspended on two 
heaters. The top view of the measurement system is shown in 
Fig. 1(f). The distance between two heaters is 10 mm. Each side of 
thermocouple is placed 1 mm from its nearest heater. Triangular 
samples with 0=20°, 40°, and 60° are studied here. 

Here we denote the positive heat flux direction when the bottom of 
rGO paper (point 1 as shown in the inset of fig. 2(a)) is at a higher 
temperature, and the negative heat fiux direction when the top (point 
2) is at a higher temperature. When the system arrives at non-equilib- 
rium stationary state, we record the input power of the heater. 
Figure 2(a) and 2(b) show the monitored temperatures as a fimction 
of time at points 1 and 2 for 60° triangular sample in the positive 
direction and the negative direction, respectively. With increasing elec- 
trical current J, the temperature systematically increase. The heat power 
caused by the electrical current through the relationship of P=rF is 
controlled to keep the same stationary state temperature difference 
between points 1 and 2, where r is the resistance of the heater. The 
most striking feature is that the magnitude of the heat power in the 
negative direction (P.) is significantly larger than that in the positive 
direction (P+). This demonstrates that there exists thermal rectification 
phenomenon in this system. Figure 2(c) shows the experimental ther- 
mal power P+ and P. as a function of temperature difference |AT| = 
I Ti — T2I . With the increase of AT, heat power monotonically increases. 
The thermal rectification coefficient R is defined as 

p 

R=^ (1) 
P+ ^ ^ 

As shown in Fig. 2(d), the rectification coefficient is 1.21 at 
|AT| = 12.7 K. The R is increased with the increasing of |AT|. The 



largest value of R is 1.28 at | AT| =39.1 K. Meanwhile, the rectification 
coefficient in rGO based asymmetric structure is much higher than that 
in carbon nanotube which is just 1.02 observed experimentally^^. The 
experimental results for 40° and 20° triangular samples can be found as 
Supplementary Fig. S4, S5 online. Figure 2(e) shows the heat power 
versus the angle of the rGO paper. The heat power P+ in the positive 
direction stays constant under different angle. Meanwhile, the heat 
power P. in the negative direction increases linearly with the angle. 
As shown in Fig. 2(f), the rectification coefficient increases linearly with 
the angle of the rGO paper. It is indicated that the rectification effect 
can be enhanced by tuning the asymmetricity. 

In order to further study the shape effect, we also measured the 
rectification of the two -rectangular shaped rGO rectifier. The dimen- 
sion of the two -rectangular shaped rGO (Wbot/Wtop= 9/2) and mea- 
surement system is shown in Fig. 3(a). Figure 3(b) and 3(c) show the 
monitored point 1 (point2) temperatures as a function of time under 
the positive (negative) temperature bias direction, respectively. We 
also found that the magnitude of the heat power in the negative 
direction (P.) is significantly larger than that in the positive direction 
(P+), which corresponds to thermal rectification effect. As shown in 
Fig. 3(d), the rectification coefficient is 1.06 at | A r| = 15.3 K. The R is 
increased with the increasing of | AT| . The largest value of is 1 . 10 at 
|AT|=46.5 K. The inset shows the experimental thermal power P+ 
and P. as a function of temperature difference | AT| . With the increase 
of AT, thermal power monotonically increases. The experimental 
results for Wbot/Wtop= 9/5 and 9/7 samples can be found as Sup- 
plementary Fig. S6, S7 online. 

Two -rectangular shaped rGO paper with Wbot/Wtop= 9/2, 9/5, 
and 9/7 are used to address the impact of width ratio. Figure 3(e) 
shows the heat power versus the Wbot/Wtop ratio of the rGO paper. 
The heat power P+ in the positive direction stays constant under 
different width ratio. However, the heat power P. in the negative 
direction increases linearly with the Wbot/Wtop. As shown in 
Fig. 3(f), the rectification coefficient increases linearly with the width 
ratio of the rGO paper. It is indicated that the rectification effect can 
be enhanced by increasing the ratio of Wbot/Wtop. 

Combining the results for both triangular and two -rectangular 
shaped graphene thermal rectifier, it is clear that the thermal rec- 
tification effect is more obvious in triangular shaped rGO paper. As 
the phonon scattering just happens at the junction of the two- 
rectangular shaped rGO, the rectification of heat power by shape 
effect is limited. Meanwhile, the phonon scattering happens along 
the whole bevel edge for the triangular shaped rGO, so it has larger 
thermal rectification coefficient. This is consistent with the idea pro- 
posed by C. Dames that interface thermal contact has negative 
impact on thermal rectification^^. 

Discussion 

To understand the underlying mechanism of the rectification phe- 
nomenon, we perform finite element simulation to study the heat 
flux and temperature distribution in the system. For the macroscopic 
asymmetric shape thermal rectifier, the Fourier's law can be written 
as: 

, dT(x) 



J=-k[xJ{x) 



dx 



(2) 



Where, x is the position, T{x) is the temperature at that position, k is 
the thermal conductivity of the rGO. 

Our simulations are based on a model that the rGO paper is in 
contact with high and low temperature heat baths. The dimensions of 
rGO paper are in the range of the experimental observation. In our 
simulations, we only consider the steady-state case. The temperature 
dependent thermal conductivity of the rGO paper can be written as^^: 

k{hj) = [3.8 X 10-^T+ 1.5 x 10-^7^ + 32 x (1 + 0.5//z)T-2]"^ (3) 
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Figure 2 | Testing results of the 60° triangular shaped solid-state thermal rectifier. Monitored temperatures (left axis) and current (right axis) as a 
function of time at points 1 and 2 of (a) positive direction and (b) negative direction, respectively. Measured (c) heat power and (d) rectification 
coefficient versus temperature difference I ATI. In Figure 2 (c) and (d), the first point is I ATI = 12.7°C (rj = 37.9°C, r2=25.2°C in positive direction); The 
second point is I ATI =21.6°C {Ti = 51.yC, T2=29TC in positive direction); The third point is I ATI =30.8°C (rj = 64.rC, r2=33.3°C in positive 
direction); The fourth point is I A Tl = 39. T C ( = 76. 1° C, 7^2= 37° C in positive direction). Measured (e) heat power and (f) rectification coefficient versus 
angle 0 under I ATI =40 K. In Figure 2(f), the first point is 0 = 2O°(rj = 73.4°C, r2=33.4°C in positive direction); The second point is 
0=4O°(r^ = 73.8°C, r2=33.8°C in positive direction); The third point is 0 = 6O°(r; = 77.5°C, T2=373°C in positive direction). 



where T is the absolute temperature and h is the thickness of rGO 
paper. 

Natural convection and radiation as the heat transfer mechanism 
between the system and the surrounding air are established. The 
surrounding environment is assumed to be stationary air at atmo- 
spheric pressure. The finite element modeling is used to perform a 
steady- state thermal analysis of the system to evaluate temperature 
profile. Ansys Multiphysics^^, a finite element analysis program, is 
employed. 

The simulation results of the temperature distribution under pos- 
itive heat flux direction and negative direction are shown in Fig. 4(a) 
and 4(b), respectively. Here, for the positive direction, temperatures 



at points 1 and 2 are set as 74 °C and 27.5°C, respectively. And the 
temperatures at points 1 and 2 are set as 27.5°C and 74°C respectively 
in negative direction. It is found that the temperature shows small 
gradient in wide rectangular shaped rGO and large gradient in nar- 
row rectangular shaped rGO, both in Fig. 4(a) and 4(b). As the 
thermal conductivity of rGO paper depends on temperature, it is 
indicated that the averaged total thermal conductivity k^y is different 
in positive and negative directions. We also show the temperature 
profile (see Supplementary Fig. S8 online) along the two-rectangular 
shaped sample. In the junction, there exhibits discontinuity in the 
temperature profile due to the sudden structural changes. The large 
temperature gradient region is a result of the presence of high 
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Figure 3 | Schematic structure and experimental results of the two-rectangular shaped solid-state thermal rectifier, (a) Schematic view of the devices. 
Monitored temperatures (left axis) and current (right axis) as a function of time at points 1 and 2 of (b) positive direction and (c) negative direction under 
Wbot/Wtop=4.5, respectively, (d) Measured rectification coefficient versus temperature difference I ATI under Wbot/Wtop=4.5. Inset is heat power versus 
temperature difference I ATI. In Figure 3 (c) and (d), the first point is I ATI = 15.3°C (rj = 37.8°C, T2=223°C in positive direction); The second point is 
\ ATI = 26.6° C {T 1 = 51. 2° Q T2= 24.6° C in positive direction) ; The third point is I A Tl = 37.8° C ( Tj = 64.3° C, 26.5° C in positive direction) ; The fourth 
point is I ATI =46.5°C ( 7^ = 74° C, r2=27.5°C in positive direction). Measured (e) heat power and (f) rectifying coefficient versus Wbot/Wtop under 
I ATI =40 K. In Figure 3(f), the first point is Wbot/Wtop=9/7 {Ti = 69.6°Q 
(Tj = 68.5°C, r2=28.5°C in positive direction); The third point is Wbot/Wtop = 



29.6° C in positive direction); The second point is Wbot/Wtop=9/5 
9/2 (ri = 66.9°C, r2=26.9°C in positive direction). 



resistance to heat flow at the junction, which translates to a reduction 
in the heat flux of the structure. The asymmetric geometry may lead 
to different radiation heat losses to the ambient. So the thermal 
radiation power in positive and negative direction should be com- 
pared with the total power. Based on the Stepan-Boltzmann's law 
(F = (jT^), the thermal radiation power can be written as: 

Pradiation = gT^ X S (4) 

Where, (7 = 5.67 x 10"^ Wm-^iC"^ is the Stepan constant, T is die 
absolute temperature and S is the surface area of the rGO paper. As the 
temperature is gradient from high to low, the temperature profile 
is simulated by finite element modeling. Then the thermal radiation 



power can be calculated by Eq. (4). As shown in Fig. 4(c), the thermal 
radiation power in positive directions is consistent with that in the 
negative directions. And both are much lower than the total power. It 
is indicated that the total power gap for positive and negative directions 
have littie relation with the thermal radiation. So the thermal radiation 
of rGO paper can be neglected. The rectification effect is mainly caused 
by thermal conduction. As shown in Fig. 4(d), the calculated rectifica- 
tion coefficient versus the temperature difference |AT| is in a good 
agreement with the experimental results under Wbot/Wtop=4.5. We 
also calculated the heat power and rectification coefficient versus 
Wbot/Wtop under |AT|=40 K. As shown in Fig. 4(e) and 4(f), die 
theoretical results are accordance with experimental results. 
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Figure 4 | Finite element modeling (FEM) simulation results of the two- rectangular shaped solid-state thermal rectifier. Simulation results of the 
temperatures distribution in (a) positive direction and (b) negative direction, respectively, (c) The total heat power (Experiments) and the thermal 
radiation power (Calculation) versus temperature difference \AT\. (d) Calculated and measured rectification coefficient versus the temperature 
difference I ATI under Wbot/Wtop=4.5. Calculated and measured (e) heat power and (f) rectification coefficient versus Wbot/Wtop under I ATI =40 K. 



As rGO paper is sensitive to humidity, we make a contrast experi- 
ment to verify the effect of humidity on our measurement. There is 
no obvious change before and after dry (see Supplementary Fig. S9 
online). It is indicated that humidity has little impact on our experi- 
ment results. 

From our simulation results, it is obvious that according to the shape 
effect, the temperature profile depends on the heat flux directions. In 
the temperature range of our experiment, the thermal conductivity of 
rGO paper decreases with temperature. For the positive heat flux dir- 
ection, the wide segment is at high temperature contact, which corre- 
sponds to high average temperature and leads to low effective thermal 
conductivity. In contrary, for the negative heat flux condition, the 
whole system is under a "lower" average temperature, which results 
in higher heat flux. Thus the temperature gradient dependent averaged 
thermal conductivity results in the thermal rectification. 



In the two-segment oxide based thermal rectifier^^'^^, its working 
mechanism is the different temperature dependence of thermal con- 
ductivity. In these two-component systems, the interface thermal 
contact has negative impact on thermal rectification^^. However, 
thermal interface material is still an emerging area. In this work, 
we propose to create macroscopic thermal rectifier based on rGO 
materials. This is a single component system, thus avoiding the issue 
of interface thermal resistance. Here the rectification ratio depends 
on the change of thermal conductivity with temperature. Figure 5(a) 
shows the temperature dependent thermal conductivity of rGO with 
different thickness, which are calculated according to the empirical 
fitting equation (3). When temperature is higher than about 200 K, 
thermal conductivity decreases with temperature, and the slope of 
the temperature dependent thermal conductivity increase with thick- 
ness, converges to a constant value when thickness is beyond 10 |im. 
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Figure 5(b) shows the thickness dependent rectification ratio. As 
expected, the rectification increases with thickness, and converges 
when thickness is larger than 10 jim, due to the impact of thickness 
on temperature dependent thermal conductivity. 

Furthermore, it is worth to address the advantage of rGO material 
by comparing with other common used materials. Figure 5(c) shows 
the temperature dependent thermal conductivity of 50 jim rGO 
paper, bulk cooper, 700 nm cooper thin fllm^^ single crystalline 
silicon layer and poly silicon^^. It is obvious that rGO paper has a 
high thermal conductivity as well as cooper. Moreover, around room 
temperature, it has more remarkable temperature dependence than 
those for cooper and single crystalline silicon film. Thus it is expected 
to result in larger thermal rectification ratio in our proposed rGO 
materials than those created by cooper and single crystalline silicon 
films. 

In summary, we demonstrated macroscopic thermal rectifier 
using rGO. The rectifier exhibits rectification which depends 
on the asymmetric shape. We analyzed the triangular and two- 
rectangular shaped rGO thermal rectifier experimentally. It is found 
that 60° triangular shaped sample has a rectification coefficient of 
1.28 and the rectification coefficient is increased with the increasing 
of the angle. In addition, the two-rectangular shaped sample has the 
rectification coefficient of 1.10 which could also be enhanced by 
increasing the width ratio. Using finite element modeling, the rec- 
tification effect for two -rectangular shaped samples was analyzed. 
The simulation results were in accordance with that of the experi- 
mental results. Although asymmetric heat flux is observed in the rGO 
system, it is worth to mention that it is not an ideal thermal rectifier, 
as heat flux in the low heat flux direction is not zero. According to 
electrical rectifier, an ideal thermal rectifier should "stop" heat flux in 
one temperature bias direction. This needs further work to reduce the 
heat flux in the low heat flux direction to make this device for better 



thermal management performance. Our results shed light on mac- 
roscopic phononic device with large-area graphene based materials, 
and may open up its important applications in thermal circuits and 
thermal management. 

Methods 

Sample preparation. Graphene oxide paper was prepared according to the Ref.^^ 
then it was reduced by acid in one hour according to the Ref.^^. The rGO detail 
preparation method is described in supplementary information. 

Characterizations. The surface morphology and chemical composition of rGO film 
were observed by SEM/EDX (DI-3100). The precise chemical composition and 
structure of rGO film was obtained by the X-ray photoelectron spectroscopy (XPS, 
Physical Electronics Inc.), X-ray diffractometer (XRD, Rigaku Inc.) and Raman 
spectroscopy using wavelength of 532 nm (HORIBA Inc.). 

Measure system setting up. First, the rGO samples were cut into triangular shape. 
Then rGO samples are suspended on two heaters. The two resistive heaters with 
R= 1245 Q is used to apply heat current Ih into the rectifier. The silver paste (Kyocera 
Chemical CT285) is used to establish good thermal contact between heaters and 
samples. The heaters are placed on the foam to reduce the heat leakage from the 
substrate. Two micro thermocouples with 50 |im dimension (KRC ST-50) are 
attached on two sides of the rectifier to monitor the temperature in the presence of the 
heat current. At last, a radiation shield is covered on the whole system to protect the 
heat fluctuation from the environment. 

Test methods. To explore the thermal rectification effect, first, the electrical current is 
applied to one of the heater. This establishes a temperature bias between the two sides. 
Here we denote the positive heat flux direction when the bottom of rGO paper (point 
1 as shown in figure 2) is at a higher temperature, and the negative heat flux direction 
when the top (point 2) is at a higher temperature. When the system arrives at non- 
equilibrium stationary state, we record the input power of the heater. In our 
measurements, we keep the same stationary state temperature difference |Ar| for 
both positive and negative temperature bias direction, and record the heat power in 
the positive direction (P+) and negative direction (P.), respectively. If the thermal 
rectification effect existed in the system, the value of P+ and P. should be different. 
The similar approach has been used by D. Sawaki et al. to study the thermal 
rectification in oxide materials 
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